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Abstract 

We apply the notion of asymptotic iteration method (AIM) to determine eigenvalues of the 
bosonic Hamiltonians that include a wide class of quantum optical models. We consider solutions 
of the Hamiltonians, which are even polynomials of the fourth order with the respect to Boson 
operators. We also demonstrate applicability of the method for obtaining eigenvalues of the simple 
Lie algebraic structures. Eigenvalues of the multi-boson Hamiltonians have been obtained by 
transforming in the form of the single boson Hamiltonian in the framework of AIM. 
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I. INTRODUCTION 



The study of the same problems from different point of view lead to the progress of the 
science and include a lot of mathematical tastes. An iteration technique [l-7] has recently 
been suggested to obtain eigenvalues of Schrodinger equation which improves both analytical 
and numerical determination of the eigenvalues and has been developed for some matrix 
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Hamiltonians, arising from the development of fast computers[8, 9]. Asymptotic iteration 
method (AIM) is very efficient to establish eigenvalues of the various quantum mechanical 
systems, because of their simplicity and low round off error. This method has been widely 
applied for determination of eigenvalues of the Schrodinger type equations. Encouraged by 
its satisfactory performance through comparisons with the other methods, we feel tempted 
to develop AIM to obtain eigenvalues of algebraic Hamiltonians. In contrast to the solution 
of the Schrodinger equation by using AIM including Coulomb, Morse, harmonic oscillator, 



etc. type potentials, the study of the algebraic Hamiltonians [ 1 OM 1 21] has not attracted much 
attention in the literature. Such Hamiltonians have been found to be useful in the study 
of electronic properties of semiconductors, quantum dots and quantum wells. It is evident 
that the formalism can also be developed for solving algebraic equations. 



The algebraic techniques have 



problems in a variety of fields 



11- 



3een proven to be useful in the description of the physical 



18] . In recent years there has been a great deal of interest 



in quantum optical models which reveal new physical phenomena described by the Hamil- 
tonians expressed as nonlinear functions of Lie algebra generators or boson and/or fermion 



operators 
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23|. Such systems have often been analyzed by using numerical methods, be- 



cause the implementation of the Lie algebraic techniques to solve those problems is not very 
efficient and most of the other analytical techniques do not yield simple analytical expres- 
sions. They require tedious calculations. In principle, if a Hamiltonian is expressed by boson 
operators, one could rely directly on the known formulae of the action of boson operators 
on a state with a defined number of particles without solving differential equations. Apart 
from the mentioned method, sometimes the Hamiltonians can be put in a simple form by 
using the transformation properties of the bosons. 

In this article, AIM is suggested and adapted to solve the bosonic Hamiltonians. We note 
that this has never been done before. As a particular case our model includes the solutions 
of the Hamiltonian of the multiphoton interactions and the Hamiltonian of the systems 
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of photons and bosons expressed in a single mode form. We briefly discuss the bosonic 
construction of the various Hamiltonians. These Hamiltonians are not only mathematically 
interesting but they also have potential interest in physics. 

The paper is organized as follows. In section 2, we briefly review the properties of 
boson and its differential realization. The procedure for solving a bosonic Hamiltonian in 
the framework of the AIM is presented in this section. Section 3 is devoted to illustrate 
determination of the eigenvalues of a bosonic Hamiltonian in the framework of the AIM. 
The bosonization of the physical Hamiltonians whose original forms are given as differential 
operators is discussed. As a practical example we illustrate the solution of the anharmonic 
oscillator and multiphoton interaction problems. In section 4, we introduce a technique to 
obtain eigenvalues of the two mode bosonic Hamiltonians by using AIM. We present the 
application of the AIM in order to obtain eigenvalues for a class of models describing two- 
mode multiphoton processes. Finally, we comment on the validity of our method and remark 
on the possible use of our method in the different fields of the physics. 

II. BASIC FORMALISM AND SOLUTION OF SINGLE BOSON HAMILTONIAN 

In this section, we illustrate solution of the single boson Hamiltonians, by modifying 
AIM. The usual differential realization of the annihilation operator a, and creation operator 
a + , are given by 

and they act on the state \n) : 

a + \n) = \Jn + 1 \n + 1) ; a \n) = y/n \n — 1} (2) 
with the commutation relation 

[a, a + ] = 1. 

A single boson Hamiltonian describing a physical system can be expressed as 

# = 5>M«> + ) 1 + E (3) 

where r ) i ■ is a constant. It is obvious that first part of then H is diagonal and exactly 
solvable. Second part of the Hamiltonian H includes non-diagonal terms and it is usually 
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solved by using various perturbation techniques. Our task is now to develop an AIM to 
obtain eigenvalues of H. We assume that action of H on the state |n) produce the following 
three term recurrence relation (or reduced to three term recurrence relation) such that 

\n + 2) = r n \n + 1) + s n \n) (4) 

where r n and s n = E — s' n are function of n. From the analogy of the AIM [l| it follows that 
(j4]) can be put a more suitable form in order to obtain eigenvalues E of H . Reformulation 
of (j3J provides the following equations: 

n = 0; |2)=r |l) + s |0)=po|l)+gb|0) 
n = 1; |3)=r 1 |2) + s 1 |l)=p 1 |l> + gi|0> 

(5) 

n = m; \m + 2) = r m \m + 1) + s m \m) = p m |1) + q m |0) 
where p m and g m are given by 

Pm T'mPm—l ^mPm—2 

Qm ^mQm—l SmQm—2 (6) 

with the initial conditions 

P-i = q~2 = 1 and p_ 2 = <Z-i = 0. 

To this end we assume that m is large enough and the states reach their asymptotic values. 
Thus we can write 

|m + 2) = p m \l)+q m \0) 

\m + 3) = p m+ i |1) +q m +i |0) (7) 
After all we can concisely write that 

Pm Pm+1 n / n 

— = or q m p m+1 - q m +iPm = 0. (8) 

qm qm+i 

The last equation can be solved for eigenvalues E, then the last approximation leads to the 
determination of the eigenvalues of the Hamiltonian H. Before going futher, we note that 
eigenvalues of the associated problem can be obtained by using the following MATHEMAT- 
ICA program code. Let us define \n) = f[n] then 



k = 20; Do[f[n+2] = Simplify [r n f[n+l]+s n f[n]], {n, 0, k}] 

(*where k is number of iteration*) 
NSolve [Coefficient [f[k + 2], f[0]]*Coefficient[f[k ], f[2]] - 

Coefficient [f[k + 2], f[2]]* Coefficient [f[k ], f[0]] == 0, El] 

(*E1 is eigenvalues of the H*) 
In the next sections, we want to illustrate our task on an explicit example. 

III. EIGENSTATE OF THE SINGLE BOSON HAMILTON! ANS 

In this section we study the determination of the single and multi-boson Hamiltonians in 
the framework of the AIM. 

1. Anharmonic oscillator 

The solution of the Schrodinger equation including anharmonic potential has attracted 
a lot of attention, arising its considerable impact on the various branches of physics as 
well as biology and chemistry. Besides its importance in physics, biology and chemistry, in 
practice anharmonic oscillator problem is always used to test the accuracy and the efficiency 
of the unperturbative methods. In this section we take a new look at the solution of the 
anharmonic oscillator problem through the modified AIM. The equation is described by the 
Hamiltonian: 



where a is a constant. Our task is now to demonstrate that the Hamiltonian (jHJ) can be 
expressed in terms of the bosons. One way to express the Hamiltonian H with boson 
operators is to use an appropriate differential realization of bosons. Using the realization 
(JT]), the Hamiltonian (jUJ) can be written as: 



H 



dx 2 



+ x 2 + ax 4 . 



(9) 




(10) 
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When the Hamiltonian (jTUjl acts on the state |n), the eigenvalue equation H\n) = E \n) can 
be transformed to the following recurrence relation: 

( H - E) \n) = (2n + 1 - E) \n) + ^-(n + n 2 + ^] \n) + 

a 



2 V 2 



V / (^ + l)(n + 2) + |n + 2) + a^/n^ - 1) - ij |n - 2} + 



-\/{n + l)(ra + 2)(n + 3)(n + 4) |n + 4) + - - l)(n - 2)(n - 3) \n - 4) = 0(11) 

Here, the skill is to express the n th even state in terms of |0) and |2) states and n th odd state 
in terms of 1 1) and |3) states. Applying the technique given in the previous section, we can 
obtain the following expressions: 

n = 0; |4> = p |0> + q |2> 
n = 2; |6> =p 2 \0) +q 2 \2) 

(12) 

n = to; \m + 4) = p m |0) + g m |2) 

n = m + 2; \m + 6) = p m+2 |0) + g m+2 |2) . 

The truncation of the state for large values of to leads to the following relations 

q m Pm+2 ~ PmQm+2 = 0. (13) 

Here Pi and qi can be calculated by using the following MATHEMATICA program code 
(again we define \n) = f[n]) 

si = Collect[Simplify[Solve[(H-E)f[n] == 0, f[n + 4]]], {%_]}] 

(*f[n+4] is obtained from (Till*) 
k = 20; Do[f[n+4] = Simplify [si [[1,1,2]]], {n, 0, k}] 

(*where k is number of iteration*) 
Solve [Coefficient [f[k + 4], f[0]]* Coefficient [f[k + 2], f[2]] - 

Coefficient [f[k + 4], f [2]]* Coefficient [f[k + 2], f[0]] == 0, El] /.«-»■ 0.1 
(*gives eigenvalues of the even state*) 
Solve [(Coefficient [f[k + 3], f[l]] Coefficient [f[k + 1], f[3]] - 

Coefficient [f[k + 3], f[3]]* Coefficient [f[k + 1], f[l]]) == 0, El]/, a 0.1 
(*gives eigenvalues of the odd states*) 
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It is obvious that, the program can easily be adapted for similar problems. The method 
introduced here gives accurate results for bosonic Hamiltonian ([TO]) . The results are given 
in Table I. As shown in the Table I our data confirm some previous results. Note that the 
results are obtained after 20 iteration. 



n 


■^present 


Em 


E [28, 29] 





1.065286 


1.065286 


1.065286 


1 


3.306872 


3.306871 


3.306872 


2 


5.747959 


5.747960 


5.747959 


3 


8.352678 


8.352642 


8.352678 


4 


11.09860 


11.09835 


11.09860 


5 


13.96993 


13.96695 


13.96993 



TABLE I: The comparison of eigenva 



numerical integration method 
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ues of anharmonic oscillator computed by the AIM [1[ , direct 



D. 



291 ] and by the present work, ATEM when a = 0.1. 



In the following subsections, it is shown that this asymptotic approach opens the way to 
the treatment of single boson quantum optical systems. 

2. A simple multiphoton interaction Hamiltonian 



Hamiltonian of t 
medium is given by 



le sing le mode coherent light with an optically bistable two photon 

24- 



26] 



H = ua + a + k [a +2 — a 2 ) + fla 



■v 



(14) 



where u is frequency, and k and Q are real constants. Time development of the Hamiltonian 
(j!4p was studied by |24| . Here we study the determination of the eigenstate of the equation 
H\n) = E \n) . The action of the Hamiltonian on the state |n) can be written as 



+ ttn(n - 1) - E) \n) + k (yjn{n - 1) \n - 2) - y/(n + l)(n + 2) \n + 2)) = 0. (15) 
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Our task is now to express n th state in terms of |0) and |1) states. 



n = 0; 


|2) 


= Po|0> 


n = 1; 


13} 


= Pi |1> 


n = 2; 


|4> 


= P2|0) 



(16) 

n = m; |m + 2) = p rn |0) 

n = m + 1; \m + 3} = p m+i |0) . 

It is obvious that eigenvalues of (CHI) can be obtained for even/odd eigenstates setting 
Pm — 0/ p m +i = 0. In this case we have used the MATHEMATICA program code given in 
SECTION II. The results are given in Table II. We have checked that the Hamiltonian (IT4")) 
can exactly be solved when £1 = 0. In this case for k — eigenvalues, E = 2n + |, and we 
have obtained the same result by using the procedure given here. 



n 




K = 0.1; O = 0.1 


K = 0.1; O = 0.5 


k = 0.5; O = 0.1 





0.5 


0.00903368 


0.00665483 


0.19828087652 


1 


2.5 


1.02298633 


1.011994512 


1.52644677404 


2 


4.5 


2.23086041 


3.010484295 


2.9397418394 


3 


6.5 


3.63572596 


6.0102248594 


4.47732150123 


4 


8.5 


5.23894189 


10.010131290 


6.1677784947 


5 


10.5 


7.04117881 


15.010086374 


8.0292960950 



TABLE II: Eigenvalues E of the Hamiltonian (|14h . for cj = 1 and various values of k and Q. 

Consequently, we have shown that AIM can be applied to the determination of the 
eigenstate of the single boson system. 

IV. EIGENSTATE OF MULTIBOSON HAMILTON! ANS 

In this section we present application of the AIM in order to obtain eigenvalues for a 
class of models describing two-mode multiphoton processes. In addition to the annihilation 
operator a, and creation operator a + , we introduce the operators b and b + in Hilbert space 
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are given by 

1 / d \ , 1/9 



•fi V * +!, J ; ^7fU + s J' (17a) 

Two boson operator, a and b, obey the usual commutation relations 

[a, b] = [a,b + ] = [b,a + ] = [a + ,b + ] = 0, [a,a + ] = [b,b + ] = 1. (18) 

Following a similar method which have been developed in the previous section, we try to 

determine the eigenvalues for a general class of two-mode multiphoton models. Hamiltonian 
of such system is given by 

H = rua + a + sub + b + n{a +s b r + b +r a s ) (19) 



In this formalism when r = s the Hamiltonian (Tl9|) satisfies the SU(2) symmetry with 



where r and k are positive integers. 

In this formalism w 
the generators 0] : 

J + = a + b, J_ = b + a, J = - (a + a - b + b) . (20) 

These are the Schwinger representation of su(2) algebra and they satisfy the commutation 
relations 

[J + ,J_]=2J [J , J±] = ±J± (21) 
The fourth generator is the total boson number operator 

N = {a + a + b + b) (22) 

which commutes with the su(2) generators. The Casimir operator of this structure is given 
by 

J = J_J + + J (J + l) = -iV(iV + 2). (23) 
If we denote the eigenvalues of the operator J by 

J = j(j + 1) (24) 

It is obvious that the irreducible representations of su{2) can be characterized by the total 
boson number N = 2j. The application of the realization (T2Q|) on a set of 2j + 1 states leads 



to the (2j + 1) -dimensional unitary irreducible representation for each j = 0, 1/2, 1, .... If the 
basis states are \j,m) {m = j,j — 1, — j), then the action of the operators on the basis 
states are given by: 

Jo\j,m) = m\j,m) 

J±\j,m) = Tm)(j±m + 1) \j, mil) (25) 
C\j,m) = j(j + 1) \ j,m) . 

An immediate practical consequence of these representation of su(2) algebra is that the 
Hamiltonian ffT9l) can easily be expressed as 

H = losN + k (J'+ JL) ■ (26) 
Eigenvalue equation H \j, m) = E \j, m) can be written as 

(2ujsj — E) \j,m) + 

■l) s (m + s - j - i)\(m + s + j)\ . . , /ri _. 

-^\j,m + s}+ (27) 



(m — j — l)!(m + j)\ 



;-l) s (-m + s - j - l)!(-m + s + j)! . 

7 : 7TT7 t j, m — s) = 

(— m — j — l)\[—m + jy. 

where N = 0,1, 2, ... .In this case the state \j, m + s) can be expressed as follows, 

m = -j; \j, -j + s) = p-j \j, -j) + q-j \j, -j - s) 
m — —j + 1; |j, —j + s + l}= p- j+ i \j, -j + 1) + q- j+ i \j, -j-s+ 1) 

m = j - 1; + s-l) = pj_i | j, j - 1) + g^-i | j, j - s - 1) 

= j; | j, j + s)= Pj-i \j, j) + qj-i \j, j -s), 

boundary condition \j, — j — s) =0. The Hamiltonian fl26l) is exactly solvable when s — 1 
and a after some straightforward treatment we can show that E = 2j + 2(n — j)/t. For the 
values s = 2 and j = 3, the the results are given in Table III. 



(28) 



V. CONCLUSION 



The basic feature of our approach is to reformulate AIM for obtaining eigenvalues of the 
bosonic Hamiltonians. Furthermore the technique given here has been used to determine 
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m 


K — J_. 

rb in ) 
1U 7 


3 


K — 95 





12 


12 


12 


±1 


| (57 ± 2V6) 


| (27 ± 2>/6) 


(9 ± 2\/6) 


±2 


J (60 ± 2\/l5) 


| (30 ± 2\/l5) 


(12 ± 2\/l5) 


±3 


i (63 ± 2^6) 


| (33 ± 2^/6) 


(15 ± 2^/6) 



TABLE III: Eigenvalues of the Hamiltonian (126j) . for cj = 1, s = 2 and j = 3. 



eigenvalues of anharmonic oscillator, multiphoton interaction problem and a class of models 
describing two-mode multiphoton processes. We have shown that AIM gives accurate results 
for eigenvalue of bosonic Hamiltonians. 

As a further work the method presented here can be developed in various directions. 
Complete spectrum of the quasi-exactly solvable problems can be obtained in the framework 
of the method presented here. Since most of the quasi-exactly solvable problems can be 
expressed in terms of generators of su(l, 1) or su(2) Lie algebra, the resulting recurrence 
relation can easily be solved by using the procedure given in this paper. The suggested 
approach can also be extended for solving boson-fermion systems. Before ending this work 
a remark is in order. This extension leads to the determination of eigenvalues of various 
Hamiltonians; Jahn- Teller Hamiltonians [19[, Rabi Hamiltonian 20], Hamiltonians of the 
Bose-Einstein condensation problems. 
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